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I.  Introduction 

Polarization  characteristics  of  underwater  light  contain  useful  additional  information  on 
inherent  optical  properties  (lOP),  concentrations  and  size  distributions  of  water  constituents 
when  compared  with  standard  reflectance  data  [1-5].  In  particular,  information  on  the  state  of 
the  water  constituents  can  be  obtained  through  analysis  of  the  spectral  and  geometrical 
angular  dependence  of  the  polarized  light  components.  In  addition,  this  analysis  can  help 
assess  visibility  in  underwater  environments,  provide  interpretation  of  ocean  lidar  signals,  etc. 
These  properties  should  also  be  taken  into  account  in  the  studies  of  atmospheric  aerosols 
above  ocean  which  employ  polarization  properties  of  atmospheric  particulates  [6]. 

Although  many  measurements  of  light  scattering  in  the  seawater  have  been  made,  the 
majority  of  them  have  not  taken  in  consideration  the  changes  that  occur  in  the  linear 
polarization  of  the  light  field.  Despite  the  importance  of  polarization  for  marine  applications, 
relatively  few  in  situ  observations  of  the  oceanic  polarization  state  of  light  have  been  carried 
out,  owing  to  a  lack  of  instrumentation  and  to  the  practical  difficulties  in  achieving  reliable 
measurements.  Watermann  [7]  was  the  first  to  measure  the  polarized  light  field  under  water 
and  showed  its  dependence  on  the  solar  zenith  angle  and  the  viewing  azimuthal  angle. 
Watermann’s  pioneering  work  was  followed  by  Ivanoff  et  aL  [8]  and  Timofeeva  [9],  whose 
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studies,  combining  both  laboratory  measurements  and  underwater  measurements  in  natural 
environments,  highlighted  the  effects  of  various  particulate  properties  (both  organic  and 
inorganic)  on  the  polarization  of  the  underwater  light  field.  In  the  1980’s  and  90’s,  a  large 
amount  of  theoretical  modeling  of  the  underwater  polarization  was  carried  out.  Much  interest 
was  focused  on  the  calculation  of  Stokes  vectors  and  Mueller  matrices  to  determine  the 
scattering  properties  of  suspended  particles  and  the  effects  of  the  air- water  interface  on  the 
underwater  polarized  light  [10,  11].  Noteworthy  exceptions  to  the  theoretical  work  are  the 
experimental  measurements  obtained  by  Voss  and  Fry  [12],  who  developed  their  own  electro¬ 
optic  light  scattering  polarimeter  to  measure  the  Mueller  matrices  of  water  samples. 

In  the  last  decade,  studies  related  to  the  underwater  light  field  have  been  gradually 
emerging  again.  This  definitely  includes  the  work  of  Cronin  and  Shashar  [13,  3]  who  were  the 
first  to  provide  hyperspectral  polarization  measurements  at  various  zenith  and  azimuthal 
viewing  angles.  Their  investigations  described  the  mechanisms  of  polarization  vision  in 
underwater  environments  focusing  on  visually  mediated  tasks  performed  by  polarization 
sensitive  animals.  However,  due  to  the  fact  that  measurements  were  mainly  taken  in  clear 
waters,  the  effect  of  suspended  particles  was  not  taken  in  consideration.  Lemer  et  aL  [14], 
following  the  work  of  Cronin  and  Shashar,  reported  polarization  measurements  not  only  in 
clear  waters  but  also  in  presence  of  a  phytoplankton  bloom  (moderately  turbid  waters). 
Reported  data  showed  high  instability  of  the  percent  polarization  in  the  presence  of  Mie 
scatterers  (as  was  the  case  of  the  phytoplankton  bloom).  The  calculated  values,  while  nicely 
reproducing  the  shape  of  the  degree  of  linear  polarization  (DOP)  in  clear  water  conditions, 
didn’t  show  the  same  instabilities.  This  is  probably  due  to  the  fact  that  the  instrument  had  to 
be  manually  adjusted  to  the  desired  azimuth  and  elevation  positions;  therefore  it  wasn’t 
suitable  for  the  rapidly  changing  underwater  environment. 

Remote  sensing  of  ocean  water  provides  information  on  suspended  particles.  Among  the 
different  types  of  suspended  matter  in  ocean  water,  phytoplankton  plays  the  primary  role  in 
global  biological  production  in  the  ocean  and,  therefore,  in  the  carbon  cycle.  Remote  sensing 
measurements  of  ocean  color  are  directly  related  to  the  water  leaving  spectral  radiance  which 
depends  on  the  absorption  and  scattering  properties  of  the  suspended  particles.  However, 
phytoplankton  cells  exhibit  only  weak  polarization  effects  (because  of  the  small  index  of 
refraction  relative  to  water)  while  inorganic  particles,  which  are  strong  backscatters, 
appreciably  affect  the  polarization  signal  [1].  Having  this  in  mind,  Chami  [15]  investigated  the 
influence  of  marine  particles  on  the  polarized  radiation  exiting  the  ocean.  Using  theoretical 
modeling  he  showed  that  an  empirical-ba.sed  inversion  approach  relying  on  the  underwater 
polarized  radiance  could  retrieve  the  concentration  of  inorganic  particles  regardless  of  the 
phytoplankton  content  in  coastal  waters.  On  this  basis,  Chami  and  McKee  [16]  also  performed 
in  situ  measurements  of  the  polarization  state  of  underwater  oceanic  radiation  with  the 
purpose  of  having  direct  estimation  of  suspended  inorganic  matter  concentration  from 
remotely  sensed  data  in  coastal  waters.  REFPOL,  the  instrument  used  in  these  measurements, 
was  a  multispectral  radiometer  with  only  four  channels  centered  at  450,  650,  850  and  1650nm 
together  with  polarizers  which  rotate  in  front  of  the  detectors,  allowing  for  successive  (not 
simultaneous)  measurements  of  radiance  values. 

In  this  paper  we  set  out  to  fill  the  gaps  in  previous  studies  to  obtain  a  comprehensive 
understanding  of  the  changes  that  occur  in  the  polarized  light  in  coastal  waters.  Unlike 
previous  multiband  instruments,  we  will  obtain  underwater  angularly  resolved  hyperspectral 
measurements  of  the  DOP  in  coastal  environments  illustrating  how  the  variability  of  the  DOP 
is  connected  to  water  constituents.  Measurements  were  also  taken  in  different  atmospheric 
conditions,  to  observe  the  effects  of  diffuse  light,  rather  than  direct  sun  light  on  the  IX)P  and  a 
comparison  between  above  and  underwater  measurements  is  explored.  Finally,  we  assess  the 
consistency  between  our  measurements  and  theoretical  polarized  radiative  transfer  models 
regarding  the  influence  of  marine  particles  on  the  polarized  signal  with  particular  emphasis  on 
the  corresponding  wavelength  and  angular  dependence. 

In  Section  2,  the  needed  definitions  and  geometrical  angles  are  given  and  the  principles  of 
the  relevant  scattering  processes  are  discussed.  In  Section  3,  the  field  experiments  will  be 
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described  and  the  operating  principles  of  the  hyperspectral  polari meter  introduced.  In  Section 
4,  measurements  of  the  DOP  will  be  presented  and  discussed  together  with  corresponding  in 
situ  standard  optical  measurements,  including  total  reflectance  spectra  (GER 
Spectroradiometer)  and  water  absorption  and  attenuation  spectra  (WET  Labs  AC-9 
instrument).  In  Section  5,  field  data  are  compared  with  simulations  obtained  using  a  Monte 
Carlo  radiative  transfer  code  of  the  atmosphere-ocean  system.  Both  radiance  and  DOP  in  the 
principal  plane  obtained  with  the  Monte  Carlo  simulations  are  presented.  Because  the  solar 
angle,  surface  wind  speed  and  detector  depths  were  recorded  along  with  the  radiance  and 
polarization  value,  to  perform  these  simulations,  only  a  few  assumptions  need  to  be  made 
about  the  scattering  characteristics  of  the  water. 

2.  Theoretical  background 

The  polarization  state  of  the  underwater  light  field  is  quantified  using  the  Stokes  vector  S  =  [I 
Q  U  V/,  where  I  represents  the  energy  flux  (Wm’^nm  *),  Q  and  U  describe  the  linearly 
polarized  comi>onent  of  this  flux,  while  its  circular  polarized  component  is  described  through 
the  Stokes  parameter  V.  Except  for  circular/elliptical  polarization  next  to  the  borders  of 
SnelPs  window  (which  is  the  circular  region  above  an  underwater  observer  with  an  aperture 
equal  to  twice  the  critical  angle,  Oc  =  48.6°),  the  underwater  light  field  is  essentially  linearly 
polarized  [8].  The  DOP  is  a  measure  of  the  percentage  of  linear  polarization  and  can  be 
expressed  as: 


DOP  = 


(1) 


/ 


Oceanic  water,  especially  in  the  coastal  areas,  contains  many  suspended  particles  (inorganic 
and  organic  particles,  which  include  living  and  nonliving  varieties  and  even  air  bubbles). 
Scattering  processes  clearly  modify  the  polarization  state  and  the  angular  features  of  the 
polarized  light  {i.e.  their  dependence  on  the  scattering  angle,  the  angle  between  the  incoming 
light  and  the  direction  of  observation),  are  strongly  related  to  the  size,  shape  and  composition 
{i.e.  refractive  index)  of  their  scattering  sources. 

In  the  real  atmosphere,  Rayleigh  scattering  by  molecules  is  the  dominant  process  affecting 
the  observed  polarization  and  the  contribution  of  the  aerosols  can  be  considered  as  a  first- 
order  correction  to  the  pure  Rayleigh  atmosphere  [17].  In  the  case  of  single  Rayleigh 
scattering  by  spherical  particles,  the  DOP  of  scattered  light  is  minimal  (actually  zero  if  the 
incoming  beam  is  1()0%  unpolarized)  at  scattering  angles  of  0°  and  180°.  In  contrast,  the 
scattering  in  the  orthogonal  plane  (90°  and  270°)  eliminates  the  in-plane  polarization 
component  so  that  the  DOP  is  maximized  with  a  value  of  1.  Multiple  scattering,  however, 
depolarizes  the  signal,  causes  a  decrease  of  the  DOP  and  the  appearance  of  neutral  points  in 
the  DOP  distribution. 

However,  in  the  ocean  environment,  scattering  by  water  and  suspended  particles  is  the 
dominant  process.  Assuming  the  shape  of  each  scattering  sources  in  the  ocean  to  be  spherical, 
and  the  composition  of  each  scattering  source  to  be  represented  by  a  single  value  of  the 
refractive  index,  the  scattering  properties  of  these  particles  can  be  computed  relatively  easily 
using  Mie  theory  in  assumption  that  the  size  distributions  and  the  refractive  indices  are  known 
[18].  In  the  case  of  Mie  scattering,  the  DOP  of  scattered  light  is  also  minimal  in  the  forward 
and  backward  directions;  however,  increasing  concentrations  of  chlorophyll  affect  the  spectral 
behavior  of  the  DOP  (for  specific  spectral  intervals,  hence  the  necessity  of  a  hyperspectral 
analysis),  and  multiple  scattering  effects,  due  to  increasing  concentrations  of  minerals,  can 
shift  the  maximum  of  the  DOP  toward  100°  scattering  angles  [1].  Figure  1  shows  the  relevant 
scattering  and  geometric  angles  we  use  in  describing  our  measurements. 
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Sun 


Fig.  I .  Geometry  of  observation.  Os^  is  the  solar  zenith  angle;  9^  is  the  detector,  or  viewing 
zenith  angle;  O^ra  tbe  scattering  angle;  <p  is  the  detector  azimuthal  angle. 


3.  Instruments  and  methods 

3.  /  The  polarimeter 

Polarization  measurements  were  taken  using  a  hyperspectral  and  multiangular  polarimeter 
developed  by  the  Optical  Remote  Sensing  Laboratory  at  the  City  College  of  New  York,  NY. 
The  instrument  eonsists  of  three  Satlantie  Hyperspectral  radianee  sensors  (reeording  intensity 
at  the  wavelengths  350-800nm,  8.5®  in- water  field  of  view)  mounted  on  a  seanning  system 
eontrolled  by  an  underwater  eleetric  stepper  motor  as  shown  in  Fig.  2.  The  heads  of  the 
sensors  are  positioned  on  the  pivot  point  of  the  motor,  in  this  way  measurements  are  always 
taken  at  the  same  depth  and  the  volume  of  water  immediately  nearby  the  heads  of  the  sensor, 
which  usually  contributes  the  majority  of  the  reeorded  signal,  was  the  same  for  all  angles.  The 
motor  rotates  the  sensors  in  a  vertical  plane  in  a  specific  angular  range  which  was  adjusted 
according  to  the  solar  zenith  angle  in  order  to  cover  the  full  0-180®  range  of  scattering  angles. 
Measurements  were  taken  every  5®.  Linear  polarizers  are  attached  in  front  of  the  sensors  with 
orientations  at  0®  (vertical),  90®  (horizontal)  and  45®. 

From  these  measurements,  the  analysis  of  total  radiance  and  DOP  is  obtained  [19]. 
Briefly,  if  7^,  I90,  Us  are  the  intensity  values  recorded  by  the  three  Satlantie  sensors,  then  the 
total  intensity  (which  corresponds  to  the  Stokes  parameter  7)  is  given  by: 


7  —  7o  +  7go  . 


(2) 


while  the  DOP  is  given  by: 


(3) 


7 


By  rotating  the  sensors  relative  to  the  nadir  direetion,  the  instrument  seanned  the  angular 
features  of  the  underwater  DOP  in  a  vertical  plane  defined  by  its  azimuth  angle  relative  to  the 
sun.  The  azimuth  angular  orientation  of  the  instrument  can  be  easily  controlled  by  hand  with 
an  accuracy  of  5®.  The  initial  azimuth  angular  position  usually  corresponds  to  the  principal 
plane,  which  is  defined  by  the  sun  and  the  zenith.  During  this  cruise  all  measurements  were 
recorded  only  in  the  principal  plane  {i.e.  (p  =  0®).  We  exploited  all  the  features  of  our 
instrumentation  {i.e.  the  possibility  of  measuring  at  different  azimuth  angles  and  at  different 
depths)  during  its  initial  testing  in  chlorophyll  rich  waters  [20].  If  measurements  are  taken  in 
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the  principal  plane  (in  which  the  Stokes  parameter  U  =  2 1 45- 1  go- 1 o  =  0),  only  the  outputs  of  the 
sensors  with  polarizers  oriented  at  0”  and  90”  are  needed  for  the  calculation  of  the  DOP  and  a 
simplified  version  of  equation  (3)  can  be  used: 

DOP  =  .  (4) 

Data  arc  acquired  through  a  customized  Labview  program  which  automatically  controls  the 
rotation  of  the  electric  stepper  motor  synchronized  with  the  data  acquisition  of  the 
hyperspectral  sensors.  The  method  of  data  collection  allows  us  to  measure  simultaneously  and 
in  real-time  the  hyperspectral  radiance  values  recorded  by  all  the  sensors  as  the  stepper  motor 
rotates.  The  total  angular  sampling  time  was  approximately  10-15  minutes  depending  on  the 
integration  time  required  to  obtain  noise-free  measurements.  Even  though  the  duration  of  the 
total  acquisition  time  should  be  decreased  to  reduce  the  strong  variability  in  the  optical 
properties  of  particles  that  can  be  found  in  highly  dynamic  coastal  ecosystems,  the  variability 
of  the  data  points  recorded  with  the  WET  Labs  AC-9  package  never  exceeded  20%,  therefore 
we  assumed  no  major  changes  in  the  water  conditions  during  the  sampling  time.  Typical 
values  of  the  integration  time  ranged  between  0.2  and  2s  per  recording,  depending  on 
scattering  angles  and  atmospheric  conditions.  Several  recordings  were  taken  for  each  angle  to 
obtain  the  averaged  value  of  the  intensity. 

3.2  Field  measurements  methodology 

Data  were  collected  during  a  recent  cruise  on  the  RA^  “Connecticut”  in  the  coastal  areas  of 
New  York  Harbor  -  Sandy  Hook,  NJ  region  (USA),  on  July  21-23  2008,  The  data  reported  in 
this  paper  were  collected  at  eight  stations.  The  coordinates  of  each  station,  the  corresponding 
solar  zenith  angles  (0sun)y  the  wind  speeds  and  the  times  of  the  day  are  given  in  Table  1 . 


Table  1.  Coordinates  and  solar  zenith  angles  of  the  sampling  stations. 


Station 

Longitude 

Latitude 

Wind  speed 

Starting  time 

1 

73‘’53’658"W 

40^6’978”N 

2\±\° 

2.5m/s 

I2X 

2 

73®53’066”W 

40®27’975”N 

57±2® 

3.8in/s 

16'" 

3 

73°52’202"W 

40'73’509”N 

36±2° 

5.5m/s 

9^ 

4 

73°54’002”W 

40°26’905”N 

22±2° 

3.2m/s 

12*" 

5 

73°47’387”W 

40°29’646”N 

56±2° 

4.8m/s 

16"' 

6 

74^8’727”W 

40‘^9’758"N 

28±2° 

6.4m/s 

To® 

7 

74°08’727"W 

40°38’023”N 

2\±r 

8.1  m/s 

12JU 

8 

75W445”W 

4r30’843”N 

46±2° 

8.0m/s 

15‘" 

The  instrument  was  lowered  from  a  winch  extending  from  the  side  of  the  ship  (approximately 
3-4m  from  the  side  of  the  ship)  so  that  shadow  effects  were  minimized.  The  principal  plane 
was  approximately  normal  to  the  ship  axes  and  the  sun  was  on  the  side  of  the  instrument.  The 
entire  assembly  was  kept  1  m  below  the  water  surface  using  four  arms  with  buoys  attached 
(Fig.  2(c)).  This  configuration  also  allowed  the  scanning  system  to  always  rotate  in  a  plane 
perpendicular  to  the  water  surface.  For  normalization  purposes,  down  welling  irradiance  was 
also  recorded  with  a  Satlantic  Hyperspectral  irradiance  sensor  positioned  on  the  front  deck  of 
the  ship.  The  wind  speed  according  to  the  ship  instrument  measurements,  as  shown  in  Table  1, 
was  in  the  range  2.5  to  8.1  m/s  and  the  ocean  surface  wave  amplitude  did  not  exceed  1.2m. 
The  sky  was  clear  blue  with  no  clouds  during  the  data  acquisition  time  for  the  first  two  days 
(Stations  1-5)  and  ovcrca.st  during  the  last  day  (Stations  6-8).  Water  optical  properties  were 
measured  at  the  same  stations  by  an  AC-9  instrument  (WET  Labs,  Inc.).  Attenuation  and 
absorption  data  are  available  for  Stations  1 , 4,  5  and  7.  Concentrations  of  mineral  particles  and 
chlorophyll  were  not  measured  directly  in  the  field,  but  they  can  be  estimated  by  analyzing 
absorption  and  attenuation  spectra  measured  with  our  AC-9  instrument.  Specifically,  the  total 
scattering  coefficient  at  555 nm,  which  is  defined  as  the  difference  between  total  attenuation 
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and  absorption  coefficients,  is  a  good  approximation  for  TSS  (Total  Suspended  Solids)  based 
on  previous  studies  [21].  Minerals  concentration  then  can  be  approximated  as  TSS 
concentration  multiplied  by  0.76  [21].  Chlorophyll  concentration,  on  the  other  hand,  can  be 
estimated  as  the  elevation  from  the  baseline  at  675  nm  on  the  absorption  spectrum  divided  by 
an  in  vivo  specific  absorption  for  chlorophyll  of  0.0146m^mg  *  [22].  Understanding  that  these 
assumptions  can  strongly  differ  for  different  regions  and  waters,  we  used  them  only  as  a  rough 
approximation  for  chlorophyll  and  minerals  concentrations. 


Fig.  2.  The  underwater  in.sirument  developed  by  the  Optical  Remote  Sensing  group  at  City 
College  of  New  York,  (a)  The  instrument  on  the  deck  of  IW  “Connecticut”,  (b)  a  detail  of  the 
Satlantic  Hyperspectral  sensors,  (c)  the  instrument  under  water. 


4.  Experimental  results 

4, 1  Characterization  of  water  compositions 

Attenuation  (c)  and  absorption  (a)  curves  with  and  without  water  are  shown  in  Fig.  3, 
minerals  and  chlorophyll  concentrations  estimated  using  the  above  approximations  are 
summarized  in  Table  2.  In  Fig.  3(c)  total  absorption  spectra  (a,^)  are  shown  which  include 
water  absorption  together  with  the  spectra  from  Fig.  3(a).  As  will  be  shown  below,  a  priori 
knowledge  of  the  total  absorption  and  attenuation  spectra  are  important  in  the  spectral  analysis 
of  the  DOP. 

We  also  recorded  total  reflectance  spectra  just  below  the  water  surface,  using  a  fiber  optic 
probe  attached  to  a  GER  spectroradiometer  (Fig.  4).  As  mentioned  earlier,  the  sky  was  clear 
with  no  clouds  at  Stations  1-5  while  at  Stations  6-8,  the  sky  was  overcast  resulting  in  mostly 
diffuse  incident  light.  Station  7  was  located  at  the  entrance  of  the  Hudson  River  and  Station  8 
was  located  further  up  the  river  which  explains  the  strong  increase  of  reflectance  (especially 
for  Station  8)  due  to  an  increase  of  mineral  scattering. 

Table  2.  Minerals  and  chlorophyll  concentrations  estimations. 


Sll 

SI4 

St5 

St7 

Minerals  (mjc/l) 

2.5 

3.2 

2.0 

3.9 

Chlorophyll  (mk/1) 

4.8 

3.9 

1.3 

2.1 

4.2  Measurements  of  particle  phase  functions 

Along  with  the  absorption  and  attenuation  coefficients,  the  volume  scattering  function  (VSF) 
is  one  of  the  fundamental  lOPs  which  govern  the  propagation  of  polarized  light  in  aquatic 
environments.  VSF  measurements  were  made  with  a  custom  device  called  the  MASCOT  [23]. 
It  uses  a  658  nm  laser  diode  source  and  17  independent  detectors  to  measure  volume 
scattering  from  10  to  170°  in  10°  increments.  Sampling  rate  for  all  detectors  is  20Hz.  A 
wedge  depolarizer  installed  in  front  of  the  source  provides  unpolarizcd  incident  irradiance. 
Independent  silicon  diode  detectors  allow  resolution  of  the  VSF  without  any  moving  parts, 
with  each  detector  gain  optimized  for  its  specific  angular  measurement.  Based  on  curvature  in 
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the  VSF  at  specific  angles,  relative  nominal  scattering  intensities,  and  the  need  to  avoid  stray 
light  contamination,  detector  field-of- views  (FOVs)  are  set  at  0.8°,  2°,  3°  and  4°  for  the  10°, 
20°,  30°  and  40°  measurements,  respectively.  Detectors  at  other  angles  have  a  FOV  of  5°. 
Calibration  and  validation  of  MASCOT  VSF  measurements  were  carried  out  with 
microspherical  beads  traceable  to  NIST  standards.  Angular  weighting  functions  were 
computed  numerically  for  each  angle  measurement.  Beam  attenuation  was  used  to  correct  for 
light  losses  along  the  optical  path  of  the  MASCOT  scattering  measurements  (all  pathlengths 
nominally  ~20cm).  Normalized  particle  phase  functions  (PPFs)  are  calculated  from  the  VSFs 
values  divided  by  the  scattering  coefficient  at  650nm  (obtained  from  the  difference  between 
attenuation  and  absorption  at  650nm  measured  with  our  AC -9  instrument).  The  PPF  for 
Station  I  is  compared  with  the  standard  Petzold’s  phase  function  [24].  Specifically,  the  three 
phase  functions  corresponding  to  each  water  type  are  averaged  to  produce  a  representative 
paniculate  phase  function,  as  described  by  Mobley  [25].  Figure  5  shows  the  comparison 
between  MASCOT  and  Petzold’s  phase  functions.  Similar  results  are  obtained  for  Stations  4, 
5  and  7,  due  to  similar  water  compositions.  The  agreement  between  MASCOT  measurements 
and  Petzold’s  functions,  justifies  the  choice  of  the  latter  to  represent  the  phase  functions  of 
particulate  matter  in  the  model  described  in  Section  5. 
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Fig.  3.  Absorption  and  attenuation  spectra  recorded  with  the  WetLabs  package. 


— St  1 


Fig.  4.  GER  total  reflectance  .spectra. 
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Station  1 


Scattering  Angle,  0^  ( ”) 


Fig.  5.  Comparison  of  MASCOT  measurements  and  standard  Petzold  functions. 


4.3  Results  of  polarization  measurements 

To  begin,  we  wish  to  verify  the  polarization  state  of  the  water  under  specialized  illumination 
conditions.  Figure  6(a)  is  an  example  (obtained  when  scattering  angle  is  0®,  which 
corresponds  to  the  sensors  looking  directly  at  the  underwater  sun  light)  where  the  DOP  is 
effectively  0  since  there’s  no  difference  in  readings  between  the  three  polarization  directions. 
In  this  case,  the  impact  of  atmospheric  panicles  and  air-water  interface  on  underwater 
polarization  was  minimal,  while  for  scattering  angles  away  from  0  and  180”  (90”,  for  example) 
the  situation  changes  dramatically  (Fig.  6(b)). 

All  recorded  DOP  values  are  based  on  measurements  of  radiances  with  different  polarizers 
in  front  of  the  sensors.  The  dependence  of  these  radiances  {i.e.  vertical,  horizontal  and  45”)  as 
well  as  of  the  total  radiance  normalized  to  the  down  welling  irradiance  as  a  function  of  the 
scattering  angle  is  shown  in  Fig  7.  The  maximum  differences  are  around  100”  which 
corresponds  to  the  maximum  of  the  DOP.  Relationships  of  these  dependencies  with  the 
volume  scattering  functions  should  be  further  investigated. 
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Fig.  6.  Spectrat  dependence  of  the  signal  recorded  by  the  Satlantic  Hyperspectral  .sensors  when 
the  scattering  angle  is  0®  (a)  and  when  it’s  90®  (b).  The  instrument  is  positioned  in  the  main 
scattering  plane,  tm  below  water. 
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station  3,  X-510nm 


Fig.  7.  Dependence  of  normalized  radiance  components  on  the  scattering  angle 

In  assessing  our  polarization  measurements  for  stability,  we  took  multiple  measurements 
when  possible.  Clearly  measurement  errors  ean  oceur  due  to  the  ehanging  ocean  environment 
and  can  be  amplified  when  DOP  values  are  small.  Figure  8  shows  the  DOP  for  all  data  points 
acquired  for  a  set  of  representative  measurements,  corresponding  to  Station  1 .  Note  the  low 
variability  of  the  experimental  data  points  for  most  of  angles  underlining  the  high  accuracy  of 
the  detection  method.  We  attribute  this  to  the  positioning  of  the  entire  assembly  on  four 
buoys,  which  kept  the  sensors  approximately  at  the  same  depth  even  in  conditions  of  high 
waves.  For  scattering  angles  inside  Snell’s  window  (Le.  between  0  and  65°  scattering  angle, 
for  Station  1),  the  largest  fluctuations  induced  by  surface  waves  were  observed. 


Station  1 ,  X-550nm 


Scattering  Angle,  0^ 

Fig.  8.  Spread  of  the  data  points  acquired  during  a  set  of  measurements.  Data  are  shown  for 
Station  I ,  X=550nm. 

To  assess  our  measurements,  we  first  note  that  the  DOP  for  ocean  waters  typically  makes 
a  bell-shaped  di.stribution  as  function  of  the  scattering  angle  with  the  maximum  around  90° 
and  going  to  zero  in  proximity  of  0  and  180°.  Typical  plots  of  the  DOP  vs.  scattering  angle, 
recorded  in  the  main  scattering  plane  at  Im  depth  are  presented  in  Fig.  9.  The  maximum  of  the 
DOP  is  lower  than  half  the  value  predicted  by  Rayleigh  theory  (i.e.  it  reaches  a  maximum 
value  of  approximately  0.4  at  410nm).  Results  for  Stations  1-5  are  very  similar  due  to  similar 
water  compositions  and  are  represented  by  Station  1  and  2  in  Fig.  9.  However,  wc  note  a 
significant  reduction  of  the  DOP  at  Station  7  and  8  (Figs.  9(c)  and  9(d)).  These  can  be  traced 
to  the  diffuse  illumination  from  clouds  as  well  as  an  increase  in  mineral  concentrations 
(Station  8,  see  reflectance  curve  in  Fig.  4).  For  reference,  downwelling  irradiance  spectra  for 
the  same  stations  are  shown  in  Fig.  9(e).  An  increased  amount  of  suspended  particles  results 
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in  increased  number  of  mulliscaltering  events,  further  depolarization  and  in  turn  lower  DOP. 
We  also  note  the  shift  of  the  maximum  of  the  DOP  towards  100°  scattering  angle.  Chami  et  al 
[1]  predicted  this  effect  and  suggested  its  use  to  allow  discrimination  between  biological  and 
non*biological  constituents  which  should  be  further  verified.  In  fact,  according  to  [3]  as  the 
real  part  of  the  refractive  index  increases  (mineral  particles)  the  absolute  maximum  of  the 
DOP  decrea.ses  in  value  and  shifts  towards  higher  scattering  angles. 
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Fig.  9.  (aHd)  Plots  of  the  DOP  vs.  scattering  angle.  The  instrument  is  located  in  the  principal 
plane  Im  below  water,  (e)  Downwelling  spectral  irradiance  recorded  at  the  same  stations. 


It’s  worth  noting  that  for  remote  sensing  purposes  only  scattering  angles  in  the  range  140- 
160°  should  be  realistically  considered.  We  recorded  above  water  measurements  for  Station  4 
and  the  data  are  presented  in  Fig.  10  in  comparison  with  underwater  measurements.  For 
angles  corresponding  to  in-water  scattering  angle  less  than  145°  (since  specular  reflection  of 
sun  light  occurs  at  -147°  in-water  scattering  angle,  for  Station  4)  abnormally  high  values  of 
DOP  were  observed,  due  to  the  contribution  of  the  polarized  reflected  light.  These  values  are 
not  related  to  the  characteristics  of  the  water  components;  hence  they  are  not  shown  in  Fig. 
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10(a).  On  the  other  hand,  for  in-water  scattering  angle  more  than  145°,  the  reflection 
coefficients  for  the  in-plane  and  out-of-plane  polarizations  are  similar  [26];  therefore  the 
observed  DOP  is  not  due  to  reflection  effects  but  to  the  polarization  properties  of  the  water¬ 
leaving  radiance.  From  this  it  is  clear  that  the  DOP  corresponding  to  the  range  140-160°  will 
not  exceed  0.2.  Taking  into  account  the  fact  that  Station  4  presented  a  relatively  low 
concentration  of  minerals  and  the  DOP  decreases  for  higher  concentration  of  minerals,  we  can 
assume  this  value  of  the  DOP  to  be  close  to  the  upper  limit  which  can  be  expected  for  remote 
sensing  measurements  in  coastal  waters.  Values  of  DOP  measured  above  and  below  water 
surface  (Figs.  10(a)  and  10(b))  for  this  range  of  angle  well  correspond  to  each  other.  The 
effect  of  absorption  is  not  noticed  here  since  the  optical  thickness  of  the  water  medium  wasn’t 
high  enough  to  decrease  the  magnitude  of  the  DOP  measured  above  water  (see  Fig.  3(b))  [27], 
when  compared  with  the  underwater  DOP.  Swell  and  wind  ruffling  of  the  water  surface  affect 
much  more  the  above  water  signal  than  the  underwater  measurements,  as  shown  by  the 
comparison  of  Figs.  10(a)  and  10(b)  curves.  The  former  contain  random  sun  glint  reflections 
not  present  underwater.  Above  water  recordings  were  performed  while  the  instrument  was 
held  by  the  winch  cable  approximately  2-3m  outboard  from  the  ship,  and  Im  above  water. 


Stalk)n  4,  Above  Water  Station  4.  Below  Water 


Fig.  10.  Plots  of  the  DOP  vs.  scattering  angle  for  Station  4.  (a)  Im  above  water,  (b)  Im  below 
water.  The  instrument  is  located  in  the  principal  plane. 


The  DOP  was  also  found  to  vary  with  wavelength.  Comparisons  between  spectral 
dependences  for  measurements  of  the  DOP  taken  at  different  stations  are  presented  in  Fig.  1 1. 
In  Fig.  1 1(a)  (which  corresponds  to  Station  1),  we  observe  a  maximum  in  the  DOP  at  lower 
wavelengths.  This  region  is  dominated  by  chlorophyll  and  CDOM  absorptions,  as  can  be  .seen 
in  Fig.  3(c)  and  in  the  normalized  absorption  spectrum,  (which  is  the  total  absorption 
normalized  to  the  maximum  total  absorption)  reported  in  the  top  part  of  Fig.  1 1(a).  In  the  top 
part  of  Fig.  1 1  (a)  the  total  absorption  spectrum  {aioi)  divided  by  the  total  attenuation  spectrum 
ic,o()  is  also  shown.  Of  course  a,ot/cto^l-a)^  where  cy,  the  single  scattering  albedo,  is  a  good 
measure  of  the  amount  of  multiple  scattering.  If  cy  decreases,  multiple  scattering  events  are 
reduced  and  the  DOP  increases.  Figure  11(a)  also  shows  that  the  DOP  reaches  maximum 
values  in  the  range  0.4-0.5  at  410  and  440nm.  On  the  other  .side  of  the  spxjctrum  {i.e.  700- 
750nm),  another  maximum  appears.  This  behavior  is  consistent  with  the  absorption  spectrum. 
After  700nm,  water  absorption  starts  increasing,  minimizing  again  elastic  scattering.  The 
relative  maximum  between  600-650nm  is  also  consistent  with  the  absorption  curve.  On  both 
sides  of  this  relative  maximum,  two  minima  occur  and  the  DOP  reaches  minimum  values 
around  0.3.  The  first  minimum  is  consistent  with  the  minimum  in  the  absorption;  absorption 
decreases  and  multiple  scattering  events  increase,  depolarizing  the  underwater  light  field.  The 
second  minimum,  however,  cannot  be  directly  related  to  the  absorption  curve.  This  dip  in  the 
DOP  is  due  to  the  chlorophyll  fluorescence  in  this  spectral  interval,  which  occurs  in  addition 
to  elastic  scattering  and  which  is  unpolarized  [4,  5].  This  hypothesis  is  confirmed  by  the  data 
recorded  at  Station  7  (Fig.  1 1(b)).  In  particular,  within  the  Hudson  River,  the  concentration  of 
chlorophyll  is  very  low,  and  the  dip  in  the  DOP  is  barely  noticeable.  In  addition,  the  effect  of 
chlorophyll  fluorescence  on  the  DOP  is  further  confirmed  by  the  comparison  of  experimental 


#106439  -  $15.00  USD  Received  1 6  Jan  2009;  revised  16  Mar  2009;  aceepled  24  Mar  2009;  published  25  Mar  2009 
(C)  2009  OSA  30  March  2009  /  Vol.  1 7,  No.  7  /  OPTICS  EXPRESS  5677 


data  with  Monte  Carlo  simulations,  which  will  be  shown  in  the  next  section.  It’s  worth 
noticing  that  the  spectral  variations  of  the  DO?  and  the  effects  of  chlorophyll  fluorescence  arc 
more  pronounced  at  scattering  angles  where  the  DOP  is  close  to  maximum  and  weaker  at 
other  angles.  In  addition,  below  400nm  and  above  700nm,  because  of  the  increasing  water 
absorption,  the  absolute  values  of  radiance  become  very  small,  leading  to  big  uncertainties  in 
the  calculation  of  the  IX)P. 
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Fig.  1 1.  Spectral  dependence  of  the  DOP:  (a)  Station  I  and  (b)  Station  7.  The  normalized  total 
absorption  spectrum  (a„ofm)  and  the  total  absorption  spectrum  divided  by  the  total  attenuation 
spectrum  (a,ot/cMt)  are  al.so  shown. 


5.  Comparison  of  modeled  and  measured  data 

5.  /  Description  of  the  radiative  transfer  model 

To  assess  the  results  beyond  simple  qualitative  experiments,  the  polari meter  measurements 
were  compared  to  the  results  of  a  Monte  Carlo  vector  radiative  transfer  model.  This  is  an 
updated  version  of  the  code  used  in  Adams  and  Kattawar  [28]  and  Adams  et  al.  [29].  It  is  a 
plane*parallel  model  of  the  coupled  atmosphere-ocean  system  and  allows  computation  for  a 
number  of  different  layers  in  both  the  atmosphere  and  ocean.  The  specific  details  can  be  found 
in  the  above  references,  and  we  will  only  give  a  brief  summary  here.  The  atmosphere  is 
modeled  as  a  two-layer  system.  The  top  layer  is  assumed  to  be  Rayleigh  scattering  only,  while 
the  bottom  layer  is  a  marine  haze  model  comprised  of  Rayleigh  plus  aerosol  scattering.  The 
aerosol  layer  above  the  ocean  surface  is  modeled  using  the  MAR-I  model  published  by  the 
International  Association  for  Meteorology  and  Atmospheric  Physics  [31].  The  sea  surface  is 
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modeled  using  the  sea-slope  statistics  of  Cox  and  Munk  [32].  In  the  water  column,  we  use  the 
dissolved  and  particulate  absorption  coefficients,  the  particulate  scattering  coefficient,  and  the 
absorption  and  scattering  coefficients  of  pure  ocean  water.  Light  scattering  in  both  the 
atmosphere  and  ocean  is  described  by  the  elastic  scattering  Mueller  matrix  for  each 
constituent.  Atmospheric  Rayleigh  scattering  is  described  using  the  Rayleigh  matrix  including 
depolarization  factors.  The  aerosol  matrix  is  found  using  the  aerosol  constituent  parameters  in 
a  Mie  code  to  produce  the  Mueller  matrix.  In  the  water,  molecular  scattering  is  again 
described  using  the  Rayleigh  matrix  and  the  depolarization  ratio  for  water.  For  paniculate 
scattering,  we  use  the  averaged  Petzold’s  particulate  volume  scattering  phase  function 
reported  by  Mobley  [25].  To  find  a  Mueller  matrix  for  particulate  scattering,  we  first  define  a 
reduced  Muller  matrix  by  dividing  each  element  by  the  (1,1)  element,  the  phase  function.  The 
particulate  Mueller  matrix  is  then  obtained  by  assuming  a  hyp)erbolic  particle  size  distribution 
with  a  slope  of  -4  and  a  constant  relative  index  of  refraction  of  ri=  1 .07+0.002/,  and  calculating 
the  reduced  matrix  using  a  Mie  code.  The  result  is  that  the  reduced  matrix  is  not  much 
different  than  a  reduced  Rayleigh  matrix,  similar  to  that  reported  in  Voss  and  Fry  [12].  In  all 
cases,  atmosphere  and  ocean,  the  phase  function  and  Mueller  matrices  are  assumed  to  be 
independent  of  wavelength.  Again,  see  reference  [29]  for  the  specific  details  and  forms  of  the 
Mueller  matrices. 

The  model  was  run  for  each  wavelength  of  the  AC-9  in-water  measurements.  The 
measured  values  of  the  dissolved  and  particulate  values  of  the  absorption  and  attenuation 
coefficients  were  used  for  the  water  column.  The  measured  wind  speed  was  used  for  the 
waveslope  statistics  (Station  I:  2.5m/s;  Station  7:  8.1  m/s).  The  largest  unknown  is  the 
atmosphere.  To  determine  the  level  of  aerosols,  we  turned  to  the  surface  irradiance 
measurements  shown  in  Fig.  9(e).  The  ratio  of  the  measured  surface  irradiance  to  the 
extraterrestrial  irradiance  was  found  for  each  wavelength.  In  the  model,  the  optical  depth  of 
the  aerosol  haze  layer  was  varied  until  the  modeled  surface  to  top  of  the  atmosphere  ratio 
matched  the  measured  one.  Once  this  aerosol  optical  depth  was  found,  the  code  was  run  for 
each  wavelength  and  the  Stokes  vector  parameters  were  found  for  a  depth  of  I  m  below  the 
surface.  These  values  were  then  divided  by  the  calculated  surface  irradiance  values  to  produce 
the  normalized  radiance  reflectance. 

5.2  Comparison  between  modeling  and  measurements. 

Figure  1 2  shows  a  comparison  between  the  measured  and  modeled  values  at  Station  1  and  7, 
at  a  wavelength  of  5 1  Onm  for  the  DOP  and  the  normalized  radiance.  Figures  1 2(a)  and  1 2(b) 
compare  the  DOP  and  radiance  at  Station  I .  The  atmospheric  optical  depth  at  this  wavelength 
was  calculated  to  be  1 .34,  and  values  for  the  other  wavelengths  ranged  from  1 .0  to  2.0.  The 
radiance  values  agree  well,  but  there  are  differences,  especially  in  the  direction  of  the  sun.  The 
model  tends  to  over  predict  the  radiance  in  this  direction,  and  then  under  predict  it  at  larger 
angles.  The  most  probable  cause  of  this  discrepancy  is  the  generally  unknown  state  of  the 
atmosphere.  By  looking  at  the  shape  of  the  modeled  radiance  curve  we  guess  that  the  incident 
skylight  was  more  diffuse  then  the  model  predicts.  The  agreement  between  the  magnitudes  of 
the  measured  and  modeled  DOP  is  very  good.  At  this  shallow  depth,  one  of  the  largest 
sources  of  uncertainty  is  the  effect  of  the  wind-blown  surface.  Though  generally  small  in 
other  azimuthal  planes,  in  the  principal  plane,  on  the  side  containing  the  direct  beam  of  the 
sun,  increased  wind  speed  will  cause  the  direct  solar  beam  to  spread  into  more  angles  and  in 
general  reduce  the  maximum  DOP.  Modeled  results  for  the  same  optical  properties  show  that 
a  smooth  sea  surface  results  in  a  maximum  DOP  approximately  5-10%  higher  than  the  wind¬ 
blown  surface  at  Station  I . 

Figures  12(c)  and  12(d)  show  the  comparisons  at  Station  7,  also  at  a  wavelength  of  51  Onm. 
This  day  was  overcast,  and  the  atmospheric  optical  depth  was  calculated  to  be  8.98,  with 
similar  values  for  other  wavelengths.  Now  it  is  unlikely  that  the  atmospheric  haze  model 
correctly  predicts  the  true  nature  of  the  atmosphere  under  these  conditions,  but  we 
hypothesize  that  the  light  field  incident  on  the  ocean  will  be  sufficiently  diffuse  that  it  doesn’t 
make  a  great  difference.  The  radiance  values  again  show  good  agreement,  though  this  time 
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the  model  tends  to  under  predict  the  radiance  at  angles  toward  the  sun.  The  DOP  again  agrees 
very  well,  though  there  is  quite  a  bit  of  statistical  noise  in  the  Monte  Carlo  simulations  due  to 
the  large  atmospheric  depth  and  very  windy  conditions. 
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Fig.  12.  Comparison  of  modeled  and  measured  data  for  510  nm.  Station  I,  (a)  DOP,  (b) 
normalized  radiance.  Station  7,  (c)  DOP,  (d)  normalized  radiance. 


The  peak  in  the  DOP  curves  near  180”  is  artifacts  of  the  statistical  nature  of  the  Monte  Carlo 
method.  This  is  a  well  known  issue  when  the  solid  angles  become  very  small  near  the  zenith 
direction,  and  there  are  insufficient  up  welling  photons  traveling  into  these  solid  angles  to 
resolve  the  DOP  sufficiently.  The  peak  observed  in  the  vStation  7  data  is  particularly 
noticeable. 

Similar  comparisons  for  a  wavelength  of  676nm  are  also  presented  in  Fig.  13.  The 
radiance  curves  match  well,  similar  to  the  510nm  wavelength  with  similar  discrepancies.  But 
now  the  DOP  match  is  poor.  The  Monte  Carlo  radiative  transfer  code  didn’t  include 
chlorophyll  fluorescence  which,  as  hypothesized  above,  would  decrea.se  the  value  of  the  DOP. 


Station  1 .  X-676nm  Station  1 ,  A.»676nm 


normalized  radiance. 
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This  probably  explains  the  significant  difference  in  the  measured  and  simulated  values. 
Figure  14  summarizes  the  spectral  comparisons  of  modeled  and  measured  DOP  for  four 
relevant  scattering  angles  for  Station  1 .  Again,  the  match  is  quite  good  except  in  the  region  of 
chlorophyll  fluorescence  as  discussed  before. 


Rg.  14.  Spectral  comparison  of  modeled  (circles)  and  measured  (solid  lines)  DOP  for  four 
relevant  scaiiering  angles,  Station  1 . 


6.  Conclusions 

Measurements  of  the  polarization  characteristics  of  coastal  waters  were  acquired  during  a 
field  experiment  in  the  coastal  areas  of  New  York  Harbor  -  Sandy  Hook,  NJ  region  (USA), 
using  a  newly  developed  hyperspectral  and  multiangular  polari meter.  An  analysis  of  the 
angular  (0-180®)  and  spectral  (400-750nm)  variations  of  the  degree  of  polarization  was 
performed.  It  was  observed  that  maximum  values  of  the  DOP,  which  occurred  approximately 
at  scattering  angle  of  100®,  did  not  exceed  0.4  for  all  stations.  In  overcast  conditions  however, 
the  light  was  still  partially  polarized  but  with  the  maximum  DOP  reduced  to  approximately 
0.2. 

Spectral  dependence  of  DOP  very  well  correlates  with  the  measured  water  absorption  and 
single  scattering  albedo  (w)  spectra:  increase  of  absorption  (decrease  of  co)  corresponds  to  the 
decrease  of  the  number  of  the  scattering  events  which  means  less  depolariz.ation.  These  effects 
are  more  pronounced  at  scattering  angles  close  to  90®  where  the  IX)P  has  its  maximal  values 
and  almost  unnoticeable  at  angles  closer  to  0  and  180®.  In  the  spectral  area  of  chlorophyll 
fluorescence  we  observed  significant  decrease  of  DOP  which  is  explained  by  the  depolarizing 
effect  of  fluorescence.  Despite  windy  conditions  at  some  stations  where  waves  were  up  to 
1.2m  high,  the  angular  profiles  of  the  DOP  and  its  maximal  values  remained  similar  to  the 
values  at  the  stations  with  the  same  water  composition  but  lower  wind  speeds. 

For  remote  sensing  applications,  when  in-water  scattering  angles  are  in  the  range  of  140- 
160®,  the  DOP  does  not  exceed  0.2.  In  addition,  we  note  that  the  measured  values  above  water 
correctly  correspond  to  the  results  of  underwater  measurements.  This  result  suggests  that, 
despite  the  effect  of  the  water  surface  on  the  polarization  of  light,  the  contribution  of  the 
underwater  polarized  light  field  is  sufficiently  significant  to  affect  the  above  water  signal. 
This  makes  the  study  of  polarization  in  the  ocean  promising  for  future  improvement  of 
retrieval  algorithms  in  complex  waters,  such  as  those  founded  in  the  coastal  areas.  In  addition, 
by  performing  these  measurements  underwater,  further  noise  sources  due  to  the  wind- 
roughened  state  can  be  eliminated. 

The  agreement  between  the  Monte  Carlo  results  and  the  experimental  data  are  also  shown 
and  clearly  demonstrates  the  success  of  radiative  transfer  simulations  applied  to  the 
transmission  and  scattering  of  light  in  an  atmosphere-ocean  system.  The  shape  of  both  the 
radiance  distribution  and  the  DOP  has  been  correctly  reproduced  for  different  atmospheric 
and  water  conditions. 
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It  is  expected  that  the  ability  of  our  polarimeter  to  provide  information  about  the 
characteristics  of  the  underwater  polarized  light  field  has  great  potential  for  application  in 
radiative  transfer  problems  in  the  earth-ocean  system,  es|>ecially  if  used  in  combination  with 
other  polarization- sensitive  instruments  recently  developed;  hypcrspectral  and  multi  angular 
data  can  be  collected  very  accurately  and  in  a  relatively  short  time,  thus  changes  both  in  the 
water  and  in  the  atmosphere  can  be  avoided. 
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